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RANDOM INTERLACEMENTS AND THE GAUSSIAN FREE 

FIELD 

By Alain-Sol Sznitman^ 
ETH Zurich 

We consider continuous time random interlacements on Z'', d> 
3, and characterize the distribution of the corresponding stationary 
random field of occupation times. When d = 3, we relate this random 
field to the two-dimensional Gaussian free field pinned at the origin 
by looking at scaled differences of occupation times of long rods by 
random interlacements at appropriately tuned levels. In the main 
asymptotic regime, a scaling factor appears in the limit, which is 
independent of the free field, and distributed as the time-marginal of 
a zero-dimensional Bessel process. For arbitrary d > 3, we also relate 
the field of occupation times at a level tending to infinity, to the 
d-dimensional Gaussian free field. 

0. Introduction. In this article we consider continuous time random in- 
terlacements on Z'^, d> 3, where each doubly infinite trajectory modulo 
time-shift in the interlacement is decorated by i.i.d. exponential variables 
with parameter 1 which specify the time spent by the trajectory at each 
step. We are interested in the random field of occupation times, that is, the 
total time spent at each site of Z*^ by the collection of trajectories with label 
at most u in the interlacement point process. 

When d = 3, we relate this stationary random field to the two-dimensional 
Gaussian free field pinned at the origin by looking at the properly scaled field 
of differences of occupation times of long rods of size N, when the level u is 
either proportional to \og N/N or much larger than log N/N. The choice of 
u proportional to log N/N corresponds to a nondegenerate probability that 
the interlacement at level u meets a given rod. In the asymptotic regime it 
brings into play an independent proportionality factor of the Gaussian free 
field, which is distributed as a certain time-marginal of a zero-dimensional 
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Bessel process. This random factor disappears from the description of the 
Umiting random field, when instead uN/ log N tends to infinity. 

For arbitrary (i > 3, we also relate the properly scaled field of differences 
of occupation times of sites by the interlacement at a level u tending to 
infinity, with the Gaussian free field on IJ^. 

Rather than discussing our results any further, we first present the model 
and refer to Section 1 for additional details. We consider the spaces 
and W of infinite and doubly infinite x (0, oo)-valued sequences, with 
d > 3, such that the Z'^-valued components form an infinite, respectively, 
doubly infinite, nearest neighbor trajectory spending finite time in any finite 
subset of Z"', and such that the (0, oo)-valued components have an infinite 
sum in the case of and infinite "forward" and "backward" sums, when 
restricted to positive and negative indices, in the case of W . 

We write X„,o"„, with n > 0, or n G Z, for the respective Z'^- and (0,oo)- 

valued canonical coordinates on and W. We denote by Px, x G Z'^, the 
law on W-^ endowed with the canonical tr-algebra, under which Xn, n > 0, 
are distributed as simple random walk starting at x, and an, n > 0, are i.i.d. 
exponential variables with parameter 1, independent from the Xn, n> 0. 
We write W* for the space W modulo time-shift, that is, W* = VF/ ~, 
where for w,w' in W, w ^ w' means that w{-) = w'{- + k) for some /c E Z. 
We denote by vr*: — )• W* the canonical map, and endow W* with the 
(T-algebra consisting of sets with an inverse image under vr* belonging to the 

canonical cr-algebra of W. 

The continuous time interlacement point process on Z'^, d > 3, is a Poisson 
point process on W x R_|_. Its intensity measure has the form u[dw*) du, 
where v is the a-finite measure on W* such that for any finite subset K 
of Z*^, the restriction of v to the subset of W* made of w* for which the 
Z'^-valued trajectory modulo time-shift enters K, is equal to vr* o Q^, the 
image of Qk under vr*, where Qk is the finite measure specified by: 

(i) Qk{Xo = x) = exix), with ex the equilibrium measure of K; 
see (1.4); 

(0.1) 

(n) when exix) > 0, conditionally on Xq = x, (X„)„>o, {X_n)n>o, 
{(7n)n£Z are independent, respectively, distributed as simple random 
walk starting at x, as simple random walk starting at x conditioned 
never to return to K and as a doubly infinite sequence of indepen- 
dent exponential variables with parameter 1. 

The existence and uniqueness of such a measure V can be shown just as in 
Section 1 of [19]. The canonical continuous time interlacement point pro- 
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cess is then constructed on a space {Q,A,¥), similar to (1.16) of [19], with 
oj = Yli>o^{w*,Ui) denoting a generic element of the set Q. We also refer 
to Remark 2.4(4) which explains how Z^, d> 3, can be replaced with a 
transient weighted graph, and continuous time random interlacements on a 
transient weighted graph are constructed. 

In the present work our main interest focuses on the collection of (con- 
tinuous) occupation times 

Lx,ui^^) ='^'^ an{wi)l{Xn{wi) = x,Ui< u} for X G Z'^,n > 
j>o nez 

(0.2) 

where u) = ''^^^{w* ,Ui) £ ^ T^*iwi) = w* for each i > 0. 

i>0 

We compute the Laplace functional of this random field and show in The- 
orem 2.3 that when V is & nonnegative function on Z"' with finite support, 
one has the identity 

L I iJ I Ei9i^x^iV{x) j 

(0.3) 

where, in the above formula, / runs over the collection of subsets of the 
support of V, gi denotes the determinant of the Green function g{-,-) re- 
stricted to / X J [see (1.1)] and c/ the sum of the coefficients of the matrix of 
cofactors of the above matrix. Both quantities are positive and [see (2.13)] 
their ratio cj/gi coincides with the capacity of /, that is, the total mass of 
the equilibrium measure ej of /. We refer to (2.26) for the extension of this 
formula to the case where Z*^ is replaced by a transient weighted graph. One 
can also consider the discrete occupation times, where an is replaced by 1 
in (0.2); however, this random field turns out to be somewhat less convenient 
to handle than (-^>x,u)xeZ'* kind of questions we investigate here; see 

Remark 2.4(5). 

The continuous time interlacement point process is related to the Poisson 
point process of Markov loops initiated in [18], which later found various 
incarnations (see, e.g.. Theorem 2.1 of [3], Sections 4 and 3 of [5], and [9], 
Chapter 9 of [8]) and was extensively analyzed in [10, 11]. Heuristically 
random interlacements correspond to a "restriction to loops going through 
infinity" of this Poisson point process; see [11], page 85. It has been shown in 
Theorem 13 of [10] (see also [11], page 61) that the field of occupation times 
of the Poisson point process of Markov loops on a finite weighted graph 
with nondegenerate killing, at a suitable choice of the level is distributed as 
half the square of the Gaussian free field on the finite graph. No such iden- 
tity holds in our context when considering a fixed level u; see Remark 2.4. 
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However, and this is the main object of this article, we present limiting pro- 
cedures which relate the field of occupation times of random interlacements 
to the Gaussian free field. 

The link with the two-dimensional Gaussian free field comes as follows. 
We look at the occupation times of long vertical rods in Z'^, by random 
interlacements at properly tuned levels. Let us incidentally mention that 
the consideration of long rods in the context of random interlacements has 
been helpful in several instances, for example, Section 3 of [15] or Section 5 
of [20]. They typically have been used as a tool in the detection of long 
*-crossings in planes, left by the trajectories of the random interlacements 
at level u, and have enabled us to quantify the rarity of such crossings when 
u is small. Here the rods in question are the subsets of Z^, 

(0.4) Jy = {x= {y, k)eZ^;l<k<N} for y G and iV > 1, 

and the corresponding Z^-stationary field of occupation times is given by 

(0.5) Cy,u=Y,^^^^^ yGZ2,n>0. 

X£jy 

We choose the levels (MAr)Ar>i and (n'^)Ar>i, so that 

log log , 

(0.6) (i) UN = a ^ witha>0, (ii) ^ =o{u]^). 

The choice in (0.6) (i) corresponds to a nondegenerate limiting probability 
exp{ — ^a} that the interlacement at level ujv does not meet any given rod 
Jy [see (4.74)] whereas the choice in (0.6) (ii) induces a vanishing limit for 
the corresponding probability. 

If we now introduce the Gaussian free field pinned at the origin, or more 
precisely [see (1-29)] a centered Gaussian field {'4>y)y^z^, with covariance 
3(a(y) + a{y') — a{y' — y)), y,y' G Z^, where a(-) is the potential kernel 
of the two-dimensional simple random walk [see (1.6)] and R an indepen- 
dent nonnegative random variable, having the law BES°(T/a, ^) of a zero- 
dimensional Bessel process at time ^ starting in -y/a at time [see (1.30)] 
we show in Theorems 4.2 and 4.9 that when N tends to infinity, 

(0.7) I , ^'"i^ I converges in distribution to the flat field with value R^ 
and that 

y,ujv .«jv \ converges in distribution to the random 



(0.8) 

field {Ripy)yez2. 
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In the case (0.6) (ii) we instead find that when N goes to infinity, 
/ \ 

(0-9) I ' 1 converges in distribution to the flat field with value 1 
and that 

— ^ — converges in distribution to {iljy)y^i2. 

There is an important connection between random interlacements and the 
structure left locally by a random walk on a large torus; see [22, 23]. In this 
light one may wonder whether some of the above results have counterparts 
in the case of a simple random walk on a large two-dimensional torus. We 
refer to Remark 4.10(1) for more on this issue. Some consequences of the 
above limit results for discrete occupation times of long rods can also be 
found in Remark 4.10(2). 

In this article, we provide yet a further link between random interlace- 
ments and the Gaussian free field, by considering the occupation times of 
random interlacements at a level u tending to infinity. If (7a;)a;eZ'' stands for 
the Gaussian free field on Z"^, d > 3, that is, the centered Gaussian field with 
covariance function El-jx^x'] = g{x,x'), x,x' G Z*^, we show in Theorem 5.1 
that when u tends to infinity, 

—Lx,u\ converges in distribution toward the flat 
(0.11) 

field with value 1 

and that 

(0.12) ( _x^u_xfi\ converges in distribution toward (7^; — 7o)a;eZ'*- 
V v2n J x&'i 

We refer to Remark 5.2 for the extension of these results to the case of ran- 
dom interlacements on a transient weighted graph and to discrete occupation 
times. 

Let us say a few words concerning proofs. We provide in Theorem 2.1 an 
expression for the characteristic function of ^^^.g^d V{x)Lx^u-, with V finitely 
supported, which shows that close to the origin it can be expressed as the 
exponential of an analytic function. This identity on the one hand leads 
to (0.3); see Theorem 2.3. On the other hand, this identity underlies the 
general line of attack, which we employ when proving the limit theorems 
corresponding to (0.7)-(0.10) and (0.11), (0.12). Namely we investigate the 
asymptotic behavior of the power series representing the above mentioned 
analytic functions. The proof of (0.8) is by far the most delicate. We analyze 



6 



A.-S. SZNITMAN 



the large N asymptotics of the power series, expressing the logarithm of 
the characteristic function of X^^g^z Wiy)^y,UN close to the origin, with 
W finitely supported on Z^, and such that '^yW{y) = 0. This asymptotic 
analysis relies on certain cancellations, which take place and enable us to 
control the coefficients of the power series. In the crucial Theorem 4.1 we 
bound these coefficients, show the asymptotic vanishing of odd coefficients 
and compute the (nonvanishing) limit of even coefficients. This theorem 
contains enough information to yield both (0.8) and (0.10); see Theorem 4.2. 
Once (0.8), (0.10) are proved, (0.7), (0.9) follow in a simpler fashion and in 
essence only require the consideration of one single rod, say Jq. The proof 
of (0.11), (0.12) in Theorem 5.1 follows a similar pattern, but is substantially 
simpler. 

Let us now describe how this article is organized. 

In Section 1 we provide additional notation and collect some results con- 
cerning potential theory, the two-dimensional free field and zero-dimensional 
Bessel processes. 

Section 2 contains the identity for the characteristic functional of the field 
of occupation times in Theorem 2.1 and the proof of formula (0.3) for the 
Laplace functional in Theorem 2.3. The extension of these results to the 
set-up of weighted graphs can be found in Remark 2.4(4). 

In Section 3 we collect estimates as preparation for the study in the next 
section of occupation times of long rods in Z^. 

Section 4 presents the limiting results (0.7)-(0.10) (see Theorems 4.2 
and 4.9) relating random interlacements in to the two-dimensional free 
field. The heart of the matter lies in Theorem 4.1, where controls over the 
relevant power series are derived. 

In Section 5 we prove (0.11), (0.12) in Theorem 5.1 and provide in Re- 
mark 5.2 the extension of these results to the case of transient weighted 
graphs, and to discrete occupation times. 

Finally let us explain our convention concerning constants. We denote 
with c,c',c,c positive constants changing from place to place. Numbered 
constants refer to the value corresponding to their first appearance in the 
text. In Sections 1, 2 and 5 constants only depend on d. In Section 3, where 
d = 3, they depend on A in (3.1), and in Section 3, where d = 3 as well, 
on A and W; see (4.3). Otherwise dependence of constants on additional 
parameters appears in the notation. 

1. Notation and some useful facts. In this section we provide some ad- 
ditional notation and recall various useful facts concerning random walks, 
discrete potential theory, the two-dimensional free field and zero-dimensional 
Bessel processes. 

We let N = {0, 1, . . .} denote the set of natural numbers. When n is a 
nonnegative real number we let [u] stand for the integer part of u. Given a 



RANDOM INTERLACEMENTS AND THE GAUSSIAN FREE FIELD 



7 



finite set A, we denote by |A| its cardinality. We write | • | for the Euclidean 
norm on M'^, d>l. For A, A' C Z"^, we denote by d{A,A') = inf{|x - x'\; 
X € A, x' G A'} the mutual distance of A and A' . When A = {x}, we write 
d{x,A') in place of d{A,A') for simplicity. We write U CC Z*^, to indicate 
that [/ is a finite subset of Z'^. Given f,g square summable functions on 
li'^ we write {f,g) = XlxeZ'* ^'^^ their scalar product. When U C Z'^, 

and / is a function on U, we routinely identify / with the function on Z"^, 
which vanishes outside U and coincides with f on U. We denote the sup- 
norm of such a function with ||/||l=o{(7), and sometimes with ||/||oo) when 
there is no ambiguity. ^ 

Given U C Z'^, we write Hu = inf{n >0;Xn£U},Hu = inf{n > 1; X„ G 
U} and Tf/ = inf{?i > 0; X„ ^ [/} for the entrance time of U, the hitting time 
of U, and the exit time from U. When p is a measure on Z*^, we denote by Pp 
the measure ^^^.g^d p{x)Px, and by the corresponding expectation. So far 
Px, X G Z*^, has only been defined when d > 3; see above (0.1). When d = 1 
or 2, this notation simply stands for the canonical law of simple random 
walk starting at x, and X„, n > 0, for the canonical process. 

When d>3, we denote by <?(•,•) the Green function 

(1.1) g{x,x') = ^Px[Xn = x'] for in Z'^. 

n>0 

It is a symmetric function, and due to translation invariance one has 

(1.2) g{x,x') = g{x' — x) = g{x — x') where g{-) = g{-,0). 
Classically one knows that g{-) < g{0), and that (see [7], page 31) 

(1.3) c'{lV\x\)-'^'^^^'> <g{x)<c{lV\x\y^'^-^^ for x G Z'^. 

When K CC Z'^, the equilibrium measure of K and the capacity of K, that 
is, the total mass of ck, are denoted by 

eK{x) = Px[Hk = oo\\k{x) for x G Z'^ and 

(1.4) 

cap(K)= Px[Hk = 'x,]. 

One can express the probability to enter K via the formula 

(1.5) Px[Hk<^]=Y, 

g{x, x')eK{x') for x G Z^. 

x'£K 

We will also consider the two-dimensional potential kernel (see (1.40), page 37 
of [7], or pages 121, 122, 148 of [17]) 

n 

(1.6) a{y)= hm V Pq [Xj = 0] - Pq [Xj = y] for y G Z^ . 

J=0 
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It is a nonnegative function on W? , which is symmetric and satisfies (cf. 
Proposition P2, page 123 of [17]) 

(1.7) lim a{y + y') — a{y') = foranyyGZ^. 

y'—^oo 

In Sections 3 and 4 [see also (0.4)] we consider long vertical rods, which are 
the subsets of defined for y (il? and iV > 1, by 

(1.8) Jy = {y}xJCZ2 where J = {1,...,A^}. 

The next lemma collects limit statements concerning the potentials of long 
rods, and in particular relates the difference of such potentials to the two- 
dimensional potential kernel. 

Lemma 1.1 {d = ^,N> l,y gZ^). 

(1.9) lim — i — y ff((0, z)) = — [with z^Z and (0, z) G /. 
^ ' N 2\ogN ^ ' 27r ' ^ ' ^ 

|z|<Af 

For X = (0, z) in Jq and y G Z^, one has 

(1.10) ^g{x,x)- ^ g{x,x") = -a{y) -bN{y,z), 
x'eJo x"eJy 

where is a nonnegative function onl? x J such that 

hN{y,z)<i)y{d{z,r)) 

(1-11) 

where lim ipy(r) = for each y G Z . 

r— >oo 

Proof. Claim (1.9) is an immediate consequence of the fact that ^ ~ 
logA^, as N goes to infinity and (cf. Theorem 1.5.4, page 31 of [7]) 

3 

(1.12) g{x) ——\x\~^ asx— )-oo. 

We now turn to the proof of (1.10), (1.11). We denote by Y_ and Z inde- 
pendent continuous time random walks on Z^ and Z with respective jump 
rates 2 and 1, starting at the respective origins of Z^ and Z. So (y,Z) is a 
continuous time random walk on Z^, starting at the origin, with jump rate 
equal to 3, and the left-hand side of (1.10) equals 

3E / l{Ys=0,Zs + z€J}ds- / l{Ys = y,Zs + z£j}ds 
Jo Jo 

(1.13) independence ^ P(p[y^ ^ _ p[y^ ^ y])P[Z, + Z e J] ds 

Jo 

= 3(/l-/2), 
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(1.14) 



h= / P[Ys = d]-P[Ys = y\ds, 
Jo 

h = / {P[Ys = 0] - P[Ys = y])P[Z, + ziJ] ds, 
Jo 



and we note that the integrand in Ii is nonnegative as a direct appHcation 
of the Chapman-Kolmogorov equation at time | and the Cauchy-Schwarz 
inequahty. If we let Yk,k >0 and Tk,k >0, (with Tq = 0), stand for the 
discrete skeleton of Y, and its successive jump times, we see that for T > 0, 

r P[Ys = 0]-P[Ys=y]ds 
Jo 

= E[{Tk+i AT-TkA T)l{Yk = 0}] 

fc>0 

(1.15) 

- E[{Tk+i AT-TkA T)l{Yk = y}] 



independence 

fc>0 



E[n+i AT-TkA T]{P[Yk = 0] - P[Yk = y]). 



Observe that Tf^^i — is an exponential variable with parameter 2, which 
is independent from T^, so that for A; > 



(1.16) ak,T = E[Tk+iAT-TkAT]=E 



Tk 



POO 

<T,2 sA{T-Tk)e-^' ds 
Jo 

decreases to zero as k tends to infinity and increases to ^ as T tends to 
oo. We set Sk = Eo<i<fc^K = 0] " = v]^ ^ > 0, so that by (1.6), 
lim^Sfc = a{y). After summation by parts in the last member of (1.15), we 
find that 

(1.17) / P[Ys = 0]-P[Ys=y]ds = y2iak,T-ak+i,T)sk. 
■^0 fc>0 

Using the observations below (1.16) we see that the left-hand side of (1-17) 
tends to |a(y) as T goes to infinity, so that 

(1.18) h = ^a(y). 

As for I2, which is nonnegative due to the remark below (1.14), we see that 
3I2 < tpy{d{z, J'^)), where we have set 

fOO ^ ^ ^ 

(1.19) i^y{r) = 3 {P[Y, = 0]-P[Ys = y])P[\Zs\>r]ds for r > 0. 

Jo 
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If is plain that ipy is a nonincreasing function, which tends to zero at infinity 
by dominated convergence. This completes the proof of Lemma 1.1. □ 

We now turn to the discussion of the two-dimensional massless Gaussian 
free field pinned at the origin. For this purpose we begin by the consideration 
of the more traditional two-dimensional massless Gaussian free field with 
Dirichlet boundary conditions outside the square Ul = [— with L > 
1; see, for instance, [2]. It is a centered Gaussian field y^y,L, y £ with 
covariance function 

(1-20) E[ipy^L^y:^L] = giiy, y) for y, y' ^I?, 

where •) stands for the Green function of the two-dimensional random 

walk killed when exiting Ul 



(1.21) gL{y^y') = E, 



'-fe>o 

Writing Hq in place of ^^{o} [see above (1.1)] it follows from the strong 
Markov property and (1.20), (1.21), that for any y G 1?, 

Vy,L - Py[Ho < TuJ(po,L is Orthogonal to (^o,L- 
Hence defining for any 7 G M, 

(1.22) ^y^L{l)=V>y,L-Py[Ho<Tu^]^o,L + Py[Ho<Tu^h, y € , 
the law of the above random field is a regular conditional probability for 
the law of {fy,L) given its value at the origin (^o,L = 7- The next lemma 
will provide two possible interpretations for the centered Gaussian field we 
consider in the sequel, in terms of the two-dimensional massless Gaussian 
free field. 



Lemma 1.2. For y,y' in 1? , one has 
(1.23) 



a{y) + a{y) - a{y' - y) = lim E[{ipy^L - ^o,L){^y',L - ^o,l)] 

L— >oo 



= lim E[<^y,L{0)^y',Lm- 

L— >oo 

Proof. By (1.20) we see that 

E[{iPy,L - ^0,L)i'^y',L - 'PO,l)] 

= 5L(y,y') +5l(0,0) -5L(0,y) -5L(0,y') 

(1.24) 

= gL{y,y) - 9L{0,y) + gL{0,0) - gL{0,y') 
- {gL{y,y) - gL{y,y'))- 

From Proposition 1.6.3, page 39 of [7], one knows that for yi,y2 in Ul, 

(1.25) gL{yi,y2)= ^ P^JXt^^ = y']a(y - 1/2) - a(yi - ^2), 
y'edUL 
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so that by (1.7) and a(0) = 0, we find that 

(1.26) lira gL{yi,yi) - gL{yi,y2) = a{yi-y2) ioiyi,y2^I?. 

Coming back to (1.24) and keeping in mind the symmetry of gii','), the 
first equality of (1-23) follows. As for the second equality, we note that 

(1.27) $y,L(0) = ify^L - VO,L + Py[Ho > TuJipo,L. 

By the strong Markov property and the symmetry of <^l(-, •) one has 

(1.26) 

Py[Ho>Tu^] = {gL{0,0) - gL{0,y))/gL{0,0) < c{y)/gL{0,0), 
and by (1.20) one finds that 
(1-28) E[^Il]=9l{0,0). 

Since lim^ (7^(0,0) = oo, it follows that the last term of (1.27) converges to 
in as L tends to infinity, and the second equality of (1-23) now follows. 
□ 

We thus introduce on some auxiliary probability space 

, s 'il)y,y ^ T? , a centered Gaussian field with covariance function 
^^•^^^ E[^yi^y,] = 3(a(y) + a{y') - a{y' - y)),y,y' G Z\ 

Up to an inessential multiplicative factor ^/S, we can thus interpret ipy^y ^ 
1? , as the field of "increments at the origin" of the two-dimensional massless 
free field, or as the two-dimensional massless free field pinned at the origin. 

The last topic of this section concerns zero-dimensional Bessel processes. 
We denote by BES''(a,r) the law at time r > of a zero-dimensional Bessel 
process starting from a > 0. If ii is a random variable with distribution 
BES'^(a,r), the Laplace transform of i?^ is given by the formula (see [14], 
page 411, or [6], page 239) 

(1.30) i^[e-^«^]=exp|-^^^| for A > 0. 

We also denote by BESQ'' (a^,r) the law of i?^; this is the distribution of a 
zero-dimensional square Bessel process at time r, starting from at time 0. 

2. Laplace functional of occupation times. In this section we obtain a 
formula for the Laplace functional of the occupation times L^^u^ which 
proves (0.3); see Theorem 2.3. As a by-product we note the absence for fixed 
u of a global factorization for the field L^^u — Lq^u, ^ S Z*^, similar to that of 
the limit law in (0.8), even through. 69^ch individuSjl VciriciblG Lx^u — ^o^u is 
distributed as the product of a time-marginal of a zero-dimensional Bessel 
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process with an independent centered Gaussian variable; see Remark 2.4(2). 
Preparatory Theorem 2.1 will be repeatedly used in the sequel and shows 
in particular that the characteristic function of a finite linear combination 
of the variables Lx^u, x (^Z,'^, is analytic in the neighborhood of the origin. 
This will play an important role in Section 4. 
We denote by G the linear operator 



(2.1) 



Gf{x)= Y 9{x,x')f{x'), xeZ", 



x'&2 



which is well defined when g{x,x')\f{x')\ < oo, and in particular when / 
vanishes outside a finite set. When V is a function on Z*^ vanishing outside 
a finite set, we write GV for the composition of G with the multiplication 
operator by V, so that GV naturally operates on L°°(Z'^) (we recall that 
II • I loo denotes the corresponding sup-norm; see the beginning of Section 1). 

Theorem 2.1. IfV has support in K dcU^, then for \\V\\oo < c{K), 

(2.2) ||Gl/||Loc^ioc < 1 
and for any u > 0, 

(2.3) eLpI V F(x)L,,Jl =exp{u(y,(/-GF)-il)}. 



Proof. Claim (2.2) is immediate. As for (2.3), we note that defining 
for X G If^, u > 0, the function on W* [see above (0.1) for notation] 

lx{w*) = Y^^n{'w)l{Xn{w) = x} for any u) E with vr* (w) = u;* , 



we have the identity 
(2.4) 



i 

for w = ^ 5(wi,u,) & ^,x £Z'^,u>0. 



The interlacement point process u is Poisson with intensity measure v{dw*) du 
under P, and hence when V is supported in K <ZC and ||l^||oo < c(-fC), 
we have 



E 



exp<^ X] ^{x)Lx,u 



:E 



exp 



vy*xiR+ 



V{x)^x{w*)1{v <u} doj{w* 
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(2.5) 



: exp< u / (e'^^s 
I Jw* 



1) dd{w* 



(0.1) 



:exp|n£'ej^ exp| 



exp<| V{X,)ds\-l 



where for s > 0, u; G we have set 

Xs{w) = Xk{w) when cro(u') H h (Tk~i{w) <s< cro{w) H h cJfc({t;) 

(by convention the term bounding s from below vanishes when A; = 0), that 
is, X_ is the natural continuous time random walk on Z*^ with jump pa- 
rameter 1 defined on W+. Thus for ||V^||oo < c(ir) we find by a classical 
calculation that 



(2.6) 



n>0 



V{Xs)ds 

VCXs,)---V{XsJdsi---dSr. 



0<si<---<s„<oo 



n>0 

= cap(ir) + 5^ 5^ eKix)[{GVri]ix), 

n>lxeK 

using Pubini's theorem and the Markov property in the last step. Since V 
vanishes outside K, it also follows from (1.5) that for n > 1, one has 

(2.7) ^e,,(x)[(Gn"l](^) = E^(^')[(GF)"-H](xO = (F,(GF)'^-H). 

X x' 

As a result we see that when ||F||oo < c(J^), 

(2.8) Ee^[e^o^ v(Xs)ds Y^^y^ (GF)"-il) = {V, {I - GVyH). 

n>l 

Inserting this identity in the last line of (2.5) completes the proof of Theo- 
rem 2.1. □ 



Remark 2.2. As a staightforward consequence of Theorem 2.1, see that 
for any finitely supported real valued function V on Z*^ and u>0, the ran- 
dom variable X^a-g^d V{x)Lx^u has a characteristic function which coincides 
in the neighborhood of the origin with the exponential of an analytic func- 
tion. So this characteristic function is analytic in the sense of Chapter 7 
of [12]. Particularly with Theorem 7.1.1, page 193 of [12], one has the iden- 
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tity 

(2.9) E 



exp<i zy^ V{x)Lx. 



X 



where in the above formula S stands for the maximal vertical strip in C 
to which the function z — )• exp{uX^„>;^ z"'{V, {GV)"'~^1)} can be an analyti- 
cally extended, and ^v,u for this extension [i.e., for z £ S, exp{z V{x)Lx^u} 
is integrable, and the equality (2.9) holds]. This fact will be very helpful and 
repeatedly used in the sequel. 

We now derive an alternative expression for the right-hand side of (2.3) 
and need some additional notation for this purpose. 

For / CC nonempty, we denote by Gj the matrix g{x,x'), x,x' G /. It 
is well known to be positive definite (see, e.g.. Lemma 3.3.6 of [13]), and we 
introduce 

(2.10) 5/ = det(G/) >0, 

where the right-hand side does not depend on the identification of I with 
{1, . . . , |/|} we use. We also set by convention gj = 1, when I = cj). Further 
we introduce 

(2.11) c/ = the sum of all coefficients of the matrix of cofactors of Gj, 

and note that c/ does not depend on the identification of / with {1, . . . , |/|} 
we employ; for instance, cj/gj coincides with the sum of all coefficients of 
the inverse matrix of Gj. The above also shows that 

(2.12) c/ > when / C C Z'' is nonempty. 

We extend the notation to the case I = (p with the convention = 0. It is 
known (see [17], page 301) that 

(2.13) cap(/) = c//g/ for all /CCZ'^. 
We are now ready for the main result of this section. 

Theorem 2.3. When V has support in K CC U^, then for \\V\\oo < 
c{K), u>^, 



(2.14) E 



where Vj *== n^g/T/ (x) (Vi = 1, by convention when I = (j)). 

In addition (2.14) holds whenever V is nonnegative and vanishes out- 
side K . 

Proof. By Theorem 2.1 we know that for ||y||oo < c{K), the left-hand 
side of (2.14) equals exp{— u(y, (/ + Gy)~"'^l)}. With no loss of generality 
we assume that > 2. We identify K with {1, . . . where n = \K\, via 
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an enumeration xi, . . . ,Xn of K. Writing vi = V{x£), we see that 

n 

(2.15) {V,{I + GVrH)= Cfc,,Vdet(I + GV), 

where G stands for the n x n matrix g(xk,Xi), 1 < k,i < n, V for the di- 
agonal matrix with coefficients vi, 1 < i < n, on the diagonal, I for the 
identity matrix and C for the matrix of cofactors of I + GV. Observe that 
for 1 <k,i <n, one has 

(2.16) Cfc,, = det((I + GV)'='^), 

where (I + GV)'^'^ stands for the n x n matrix, where the A;th line and the 
ith column of 1 + GV have been replaced by 0, except for the coefficient at 
their intersection, which is replaced by 1. 

Given an n x n matrix A = (a^^^), we develop the determinant of A ac- 
cording to the classical formula 

n 

(2.17) detA = ^sign(f7) JJafc,^(fc), 

cr k=l 

where a runs over the permutations of {1, . . . and sign((T) denotes the 
signature of a. 

We now develop the determinant det(I + GV). For each subset J C 
{1, . . . , n} we collect the terms corresponding to permutations a of {1, . . . , n} 
such that a{k) = k, for k £ J, the choice of 1 in each term {l + g{0)vk), k £ J, 
and for any k ^ J such that a^k) = k, the choice of g{0)vk instead. Thus 
for each such J, setting J = {1, . . . ,n} \ J, the sum of these terms equals 
det(G| j^j) n^gj"^^- Thus summing over all subsets J of {1, ... ,n} we find 

(2.18) det(I + GV) = ^ gjVj. 

ICK 

We now turn to the numerator of the right-hand side of (2.15). We use the 
convention {k,i} = {k} = {£}, when k = L As above we develop the deter- 
minant det((I + GV)'^'^); see (2.16), (2.17). We can assume that the permu- 
tations a of {1, . . . ,n} entering the development satisfy <j{k) = I. For each 
J C {1, . . . ,n} \ {/c,^}, we collect the terms corresponding to permutations 
a such that a{rri) = m, for m £ J, the choice of 1 in each term (1 -|- g{0)vm), 
for m£ J, and for any m ^ J U {k, i} with a{m) = m, the choice of g{0)vm 
instead. Setting J = {1, . . . , n} \ J, we see that the sum of these terms for a 
fixed given J as above equals 11 r„iUmdet(G~^~), where G ~ stands 
for the matrix where the A:th line and the ith column of the matrix G^j^j 
(i.e., G restricted to J x J) are replaced by zero except for the coefficient 
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at their intersection, which is replaced by 1. Thus summing over all possible 
J C {1, . . . ,n} \ and all in {1, . . . ,n}, we obtain 

n n 

Ck,eve = Y J2 H ^™det(G^^^^) 

k,e=l k,i=lHD{k,e}m£H 

(2.19) = n E det(G^t/.) 

(l>7^HC{l,...,n}m£H k/£H 



(2A1) 

ICK 



(using the convention c,^ = in the last equality) . 

Combining (2.15), (2.16), (2.18) we obtain (2.14). Finally in the case of 
a nonnegative V with support in K, we note that ¥,[exp{—z'^^V{x)Lx^u}] 
is analytic in the strip Rez > 0, and coincides for small positive z with 
the function exp{— u ^^^^ c/V/z'^l/ ^^^^^ gf/V/zl^l}, which is analytic in 
the neighborhood of the positive half-line. Both functions thus coincide for 
z = 1, and our last claim follows. □ 

Remark 2.4. (1) Choosing V = XIk, with A > and K CC Z*^, we 
deduce from (2.14) by letting A tend to infinity that 

¥[L^ u = 0, for ah x£K] = expi -u— \ for n > 0. 

I 9K } 

Introducing the interlacement at level u, 

X^{uj) = {x G U^] for some i > such that Uj < u, w* enters x} 

i>0 

and taking (2.13) into account, we recover the well-known formula (see (2.16) 
of [19]) 

(2.20) F[VnK = (/)] = exp{-ucap{K)} for u>0,K CCZ'^. 

(2) Choosing V = Alja-}, with A > and x G Z*^, Theorem 2.3 now yields 
that 

(2.21) £[exp{-AL^,„}] = expi ^^^^rv I ^ ^ 0' 



1 + 5(0)AJ 

and in view of (1.30) we find that 

(2.22) L^^u is BESQ° (u, ^ ) distributed. 
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If x,x' G are distinct, choosing V = z{l{x} ~ ^{x'}) (2-14) with z small 
and real and extending the identity to z = it, t G M, with the help of (2.9), 
we find that 

(2 23) 

^\ l + (c/(0)2-5(x'-x)2)t2 

In view of (1.30) we thus find that 

Lx,u — Lx'^u has the law of Rip, where R and are independent, 
(2.24) respectively, BES^(-v/n; ^ ), and centered Gaussian with 

variance 4(g(0) — g{x — x')) distributed. 

Let us, however, point out that in the case of three distinct points x, x' , x" in 
Z'^, the law of the random vector {Lx'^u ~ Lx,u, Lx"^u ~ Lx,u) does not coincide 
with that of the scalar multiplication of a two-dimensional Gaussian vector 
by an independent BES^(a, r)-variable, when u > 0. Indeed one has 

¥[Lx>,u - Lx,u > 0, L^/> - Lx,u = 0] > P[X" 3 2;',X" n {x, x"} = 0] > 

as a consequence of (2.20) and the fact that cap({x, x', x"}) > cap({x,x"}). 
But for the above mentioned distribution both components necessarily van- 
ish simultaneously on a set of full measure, and the above probability would 
equal zero if such an identity in law was to hold. We will, however, see in 
Section 5 how x G Z'^, can be related to the d-dimensional Gaussian 

free field, by letting u tend to infinity, instead of keeping u fixed. 

(3) Random interlacements can be related to the Poissonian gas of Markov 
loops; see [10, 11]. Heuristically they correspond to "loops passing through 
infinity;" see [11], page 85. The identity for Markov loops corresponding 
to (2.3) of Theorem 2.1 above can be found in Corollary 1 of Chapter 4, 
Section 1 and Proposition 7 of Chapter 2, Section 4 of [11]. The presence of 
a logarithm and a trace in the expressions leading to Proposition 7 of [11] 
is emblematic of the Markov loop measure and can be contrasted with the 
expression in (2.8) for random interlacements [which is then inserted in the 
last line of (2.5)]. In the case of a Poissonian gas of Markov loops on a finite 
weighted graph with a suitable killing, it is shown in Theorem 13 of [10] that 
the occupation field of the gas of loops at level ^ (playing the role of u in the 
context of [10]) is distributed as half the square of a centered Gaussian free 
field with covariance the corresponding Green density. For similar reasons 
as in (2) above, no such identity holds for random interlacements at any 
fixed level u. We will, however, present in the next two sections limiting 
procedures that relate random interlacements to the Gaussian free field. 

(4) As we now explain, the results of this section can be extended to 
the case of continuous time random interlacements on a transient weighted 



for t G M. 
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graph. One considers a countable connected graph E which is locally finite 
and endowed with nonnegative symmetric weights px^x' = Px',x, which are 
positive exactly when {x, x'} belongs to the edge set £ of E. One assumes 
that the induced random walk with transition probability Px,x' = Px,x' / pi^), 
where p{x) = Ylx'^E Px,x' ■, is transient. Random interlacements can be con- 
structed on such a transient weighted graph; see [19], Remark 1.4, and [21]. 
Continuous time random interlacements can also be constructed, in essence 
by the same procedure described in the Introduction, endowing the discrete 
doubly infinite paths with i.i.d. exponenticd variables of parameter 1. The 
corresponding expression for the measure Qk^ for K finite subset of E, re- 
mains the same as in (0.1); simply, the expression for exi') the equilibrium 
measure of K, which appears in (1.4), now has to be multiplied by the factor 
p{x) in the present context. 

The occupation time variables Lx^m x G E, u>0, are defined by a similar 
formula as in (0.2), but the expression on the right-hand side of (0.2) is now 
divided by p{x). The linear operator G corresponding to (2.1) operates, say, 
on functions f on E with finite support, via the formula 



Of(x) 



E 

x'eE 



g{x,x')fix')pix'), xeE, 



where g{-, •) now stands for the Green density, which is obtained by dividing 
the expression corresponding to the right-hand side of (1.1) by p{x'). 

The proof of Theorem 2.1 can be adapted to this context to show that 
when K is a finite subset of E, and V has support in K, then for ||V^||loo(^) 



sufficiently small, ||Gy||io 



(2.25) 



E 



exp<^ X] ^ix)Lx,uP{x) 



{E)-^L°°{E) < li and for any n > 0, 

■.exp{u{V,{I-GV)-h)}, 



x£E 



where now {f,g) stands for Xlxs-B (whenever this sum is abso- 
lutely convergent). Likewise the proof of Theorem 2.3 is easily adapted, and 
one finds that for V as above, n > 0, 

/ \T 11 f E/cxC/n:,.e/l/(x) 

eyip<-2_^V{x)Lx,u} 

^ x&E ^ - 



(2.26) E 



■ exp 



T,ick9i'^x&iV{x) 



with gj and cj defined as in (2.10), (2.11) [with g{-,-) now denoting the 
Green density]. 

(5) One can define the stationary field of discrete occupation times 
x G Z*^, u>0, analogously to Lx^u, simply replacing ct„ by 1 in (0.2). When 
F is a function on Z'^ with support contained in K CC Z'^, it follows that 
1 — is a function supported in K with values in (— oo, 1), and one has 
the identity 

(2.27) E exp|^y(x)4,«l =E exp 



5:(l-e-^(^))L,„ 



u>0, 
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as can be seen by integrating out the exponential variables in the right 
member of (2.27) (of course both members of the above equality may be 
infinite). As a result, Theorems 2.1 and 2.3 also yield identities concerning 
the Laplace functional of i^x,u)xezd-- 

3. Preparation for the study of long rods. In this section we introduce 
notation specific to Z'^ and provide estimates in Lemmas 3.1 and 3.2, which 
will be recurrently used in the next section, when we investigate the occu- 
pation times spent by interlacements at a suitably scaled level in long rods. 
These controls will play an important role in the asymptotic analysis of the 
power series entering the characteristic functions of these occupation times. 
Throughout this section we assume that d = 3, and constants depend on the 
finite subset A of introduced in (3.1) below. The notation || • ||oo refers to 
the supremum norm || • where B appears in (3.1). 

We consider A C C containing and > 1 . We also define 

(3.1) B = AxJ where J = {1,...,A^}. 

We write tt^^ and ttz for the respective Z^- and Z-projections on Z^ identified 
with Z^ X Z. Given a function F on B, we write (F) for the function obtained 
by averaging F on horizontal layers and {F)z for the average of F on the 
layer A x {z}, so that 

(3.2) (F)(x) = (F), = ^ J];F((y,z)) for x G 5 with 7rz(x) = 
We also introduce the function 

(3.3) [F]oix) = F((0, z)) for xinB with 7rz(x) = z. 
It is plain that for any function F on B, 

(3.4) {F-{F)) = 
and that 

(3.5) F = (F) when F only depends on the Z-component. 
We consider nonempty sub-intervals of J, 

(3.6) h^h^J with L = d(/o, J\/i) > 1 
(we refer to the beginning of Section 1 for notation). We write 

(3.7) Co = A X Jo C = A X /i C 

We recall the convention concerning constants and the notation || • ||oo stated 
at the beginning of this section. The estimates in the next lemma reflect the 
decay at infinity of the Green function [see (1.12)] and the fact that the 
discrete gradient of g{-) has an improved decay at infinity; see (3.15) below. 
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Lemma 3.1. For any function F on B, one has: 

(3.8) ||G'F||oo<clogiV||F||oo; 

(3.9) \\Ic,GIb\c,FU < clog(^^) \\F\U, 

(3.10) IIGFIU <c||F|U when{F) = 0; 

(3.11) \\IcoG1b\c,F\\^<^\\Ib\c,F\\oo when{F)=0; 

(3.12) \\GF-[GF]o\\^<c\\F\U, 

(3.13) \\1co{G1b\c,F -[GIb\c,F]o)\\oo < ^\\Ib\c,F\\^- 

Proof. We begin with (3.9) and note that for x £ Co, 



\{GIb\c,F){x)\ 



g{x,x')F{x') ^<\( l)\\F\\oo 

x'eB\Ci ^L<k<N ' 

<clog(^)||F|U, 

whence (3.9). The bound (3.8) is proved in the same fashion. 

We then turn to the proof of (3.11) and note that when {F) = 0, for x G Co 
one has with the notation x' = {y', z'),x = {y, z') [so that ttz{x') = irzix) = z'\ 

(G\b\c,F){x) = Y 9{x,x')F{x')= Y E^(^'^')^(^') 
(3.14) <a=° 1 ^ <7(x,x')(F(x')-F(x)) 

' ' z'eJ\hy',yeA 
From Theorem 1.5.5, page 32 of [7], one knows that 



c a 



(3.15) \g{x + a) - g{x)\ < Y^^-y-^ for x G Z'^, |a| < diam(A), 
where diam(A) stands for the diameter of A. As a result we see that 

(3.16) \iGlB\c^F)ix)\ < c\\1b\c,F\\^ ^ TT^ " zH^^Vt^i^H^' 

k>L 

whence (3.11). One proves (3.10) analogously. 
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As for (3.13) we note that with the notation xq = (0,z) when x = {y,z), 
we have for x G Co, 



\{G1b\c-,F-[G1b\c-,F]o){x) 



{gix, x') - g{xo,x'))F{x') 

x'eB\Ci 



(3.15) 1 



L>k 
C 

whence (3.13). The bound (3.12) is proved analogously. □ 
We conclude this section with the following lemma. 

Lemma 3.2. For F,H functions on B, one has 

(3.18) \\{F{GH))\\^<c\\F\\^\\H\\^ when {F) = or {H) = 0. 

Proof. The case (H) = is immediate thanks to (3.10). In the case 
where (F) = 0, we write for z £ J, with the notation x = {y,z),x = (jj, z), 

{F{GH)). = ±S2n^)T.9ix,x')H{x') 

' ' j/eA x'&B 

(3.19) ^^-'W ^ iF{x)-F{x))Y,9{x,x')H{x') 

' ' y,j/6A x'eB 



^ {g{x,x')-gix,x'))F{x)H{x'). 

' ' y,y&A,x'&B 



By (3.15) we thus find that with the notation 7rz{x') = z' , 

y,yeA,x'&B 

(3.20) 



< r F\\ H\\ 
and (3.18) follows. □ 

4. Occupation times of long rods in Z^. In this section we relate the field 
of occupation times of long rods in by random interlacements at a suitably 
scaled level [see (4.1) below] with the two-dimensional free field pinned at 
the origin introduced in (1.29). The main results are stated in Theorems 4.2 
and 4.9. The approach is roughly the following. By Theorem 2.1 we can 
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express the characteristic functionals of the scaled fields of occupation times 
of the long rods as exponentials of certain power series. The main task 
is to control the asymptotic behavior of these power series. This analysis is 
carried out in the central Theorem 4.1 as well as in the simpler Theorem 4.8. 
Throughout this section we assume that d = 3. The constants depend on the 
finite subset A of [cf. (3.1) and above (4.3)] as well as on the function W 
with support in A that appears in (4.3). As in Section 3 we denote by || • ||oo 
the supremum norm || • ||loo(b), with B as in (3.1). 

We consider a > 0, and a positive sequence 77V tending to infinity. We will 
analyze the random fields of occupation times of the long rods Jy, y € Z^, 
with Jy = {y} X J = {y} x {1, . . . , N} C Z^, by random interlacements at the 
scaled levels 

log N log N 

(4.1) UN = a ^ and Uj^ = "fN ^ withiV>l. 

The corresponding occupation times of the rods Jy, y G Z^, are 

Let us point out that sequences of levels converging faster to zero than un 
are not interesting in the present context; see Remark 4.3 below. 

As in Section 3, we consider some A CC Z^ containing 0. Further we 
introduce a function W on 1? such that: 

(4.3) (i) VF(y) = outside A; (ii) ^W{y) = Q. 

y 

We define the functions on Z^, 



(4.4) 



so that 



Vn{x) = — l=VF(y)lj(^), 
VlogiV 

V'^{x) = -^Vn{x) withx = (y,z), 



(4.5) C^"^ ^(y)^^ = E yN{x)L^,u. 



y&? 



and similarly 



(4.6) C!^ ^(^) A% = E ^i^(^)^^,nV 
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It follows from Theorem 2.1 and Remark 2.2 that 

(4.7) E[exp{2;£Ar}] =exp< ^^ajv(n)z" > for \z\ < r^r in C 

with r/vT > and 

(4.8) aN{n)=UN{VN,{GVNT~^l) forn>l. 
As a result of the centering condition (4.3) (ii) we have 

(4.9) ajv(l)=0. 
In a similar fashion we have 

(4.10) E[exp{z£'^}] = expl J]] a'^(n)z" I for \z\ < in C 
with r'j^ > and 

(4.11) aM = ^M, (Gy^)"-H) ^'-'U'-'^ -7^""^'«Jv(n) for n > 1. 

a 

The heart of the matter for the proof of Theorem 4.2 lies in the analysis of 
the large N behavior of the coefficients ai\i{n), n > 1. There is a dichotomy 
between the case of odd n, with an asymptotic vanishing of aN{n), and 
even n, with a positive limit of aAr(n), as N goes to infinity. The crucial 
controls are contained in the next theorem. We recall the convention con- 
cerning constants stated at the beginning of this section. 

Theorem 4.1. 

(4.12) \aN{n)\<ac^ for all n>l,N > 1, 

(4.13) foranyk>0 limaAr(2A; + 1) = 0, 

where we have set [see (1.6) for notation] 

(4.15) £{W) = -3Y,W{y)W{y')a{y' -y). 

y,y' 

Note that due to (4.3) (ii) we can express £{W) in terms of the two- 
dimensional Gaussian free field tpy, y G Z^, introduced in (1.29), via the 
formula 

(4.16) £iW) = E\(Y,W{y)iPy^ ' 

Before turning to the proof of Theorem 4.1, we first explain how this theorem 
enables us to derive the convergence in law of the appropriately scaled fields 
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Cy^N — ^o,N, y G ^ and C'^ — C'q j^, y G Z . We tacitly endow M with 
the product topology, so that the convergence stated in Theorem 4.2 actu- 
ally corresponds to the convergence in distribution of all finite-dimensional 
marginals of the relevant random fields. The main result of this section is 
the next theorem, which proves (0.8) and (0.10). 

Theorem 4.2. 

(4 17) ^'^ ^ ^^^^ infinity, ( ^.yiogjv' ~-^?^£^^' '^'^^^W^^ distribution 
to the random field {Ripy)y^^2 , 

where R and {ipy)y£'i'^ are independent and 

(4.18) R is BES°(Va, distributed, 

(4.19) i4'y)y(£Z^ is the centered Gaussian field introduced in (1.29). 
Moreover, 

(4 20) '^^ ^ ^"^^ infinity, ( ^■J——^ —)y^^2 converges in distribution 

to ('0S/)yeZ2- 

Proof (assuming Theorem 4.1). We begin with the proof of (4.17). We 
consider W as in (4.3) and Cn as in (4.5). By (4.7), (4.12) and Remark 2.2, 
we see that exp{zCN} is integrable for any z with cq\ Rez| < 1, and that 

(4.21) E[exp{z£7v}] =exp< ^^aAr(n)z"' > for any \z\ < Cq^ in C. 
In particular (4.12) implies that 

(4.22) supE[cosh(r£7v)] < oo when r < Cq ^. 

TV 

As a result the laws of the variables C]\f are tight, and the variables exp{z£Ar}, 
> 1, with |Rez| < r < Cq ^, are uniformly integrable. If along some subse- 
quence Nk, k>l, the variables Cnj^ converge in distribution to £, it follows 
from Theorem 5.4, page 32 in [1], that for \z\ < Cq ^, 

£'[exp{z£}] = limE[exp{2£Arj.}] 

(4.23) =limexpj^a^,(n)z"| 

= exp!.-£{W)z^/(l-—£{W)z'^ 



^2 ' ' / \ 27r 

using Theorem 4.1 in the last equality. This determines the characteristic 
function of the law of C, and by (1.30) shows that S has same distribution as 
Rip where i?,^ are independent variables with R'BES^{^/a, ^)-distributed. 
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and ip a centered Gaussian variable with zero mean and variance £{W). This 
proves that for any W as in (4.3), Cn converges in distribution to Rip as 
above, when N tends to infinity. In view of (4.16), this completes the proof 
of (4.17). 

The proof of (4.20) is analogous. Due to (4.11), we know that a'j^{n) = 
«^iv "''^OAr(n), and in particular a'jy(2) = ^aAr(2) converges to ^£{W), 
whereas for n 7^ 2, a'j^{n) converges to zero as N goes to infinity. We can 
use the same arguments as above and find that C'j^ converges in distribution 
to a variable C such that for small |z| in C, E[exp{z/:'}] = exp{i^(VF)z2}. 
The claim (4.20) then follows immediately. □ 

Remark 4.3. For the kind of limit theorems discussed here, sequences 
of levels converging to zero faster than um lead to trivial results, as we now 
explain. In a standard way (see, e.g.. Remark 3.1(3) in [16]), one has the 
bound cap(Jj^) < c ^^^^ , for all y G 1? . If we pick u'^ so that u'^ = o( '°|/^ ), 
then (2.20) implies that for all y G Z^, with probability tending to 1 as 
N goes to infinity, the interlacement at level u"^ does not intersect Jy. In 
particular, if we define Cy ^ in analogy to Cy^N in (4-2) with u"^ in place of 
UN, the random field {^y N)y^z'^ converges in distribution to the constant 
field equal to zero, as goes to infinity. 

Proof of Theorem 4.1. We recall the convention concerning con- 
stants and the notation || • ||oo stated at the beginning of this section. The lin- 
ear operators under consideration throughout the proof will be restricted to 
the space of functions vanishing outside B. The centering condition (4.3) (ii) 
and the ensuing identity (Vn) = play a crucial role. We will first prove (4.12), 
(4.13). Our first step is to control the norm of the operators (GVn)'^- 

Lemma 4.4. 

(4.24) \\{GVNf\\L^(^B)^LooiB)<Cl. 

Proof. Given a function F on B, we write in the notation of (3.2), (3.3) 
{GVN)^F = Ai + A2 + A-i where 
Ai = GVnG{VnF-{VnF)), 

(4.25) 

A2 = GVn[G{VnF)]o, 
As = GVn{G{VnF) - [G{VnF)]o). 
With the help of Lemma 3.1 we see that 

(3.8), (4.4) (3.10), (4.4) 

(4.26) pilloo < cy^h^\\G{VNF-{VNF))\\^ < c\\F\\^. 
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Since {V]\f[G{VNF)]Q) = 0, we also find that 

(3.10), (4.4) c (3.8), (4.4) 

(4.27) P2II00 < -^==\\G{VnF)\\^ < c\\F\ 
Finally we have 

(3.8), (4.4) 



P3II00 < c^/logN\\G{VNF-[G{VNF)]o)\\^ 
(4.28) 

(3.12), (4.4) 

< c||F||oo. 

Collecting (4.25)-(4.28), the claim (4.24) now follows. □ 

Before proving (4.12), (4.13) we still need the following lemma, which 
shows that the kernel of the linear operator F — )• (VnF) , is almost invari- 
ant under (GVn)'^. We will later see [cf. (4.61)] that the function F = 1b, 
which belongs to this kernel, is in an appropriate sense, close to being an 
eigenvector of (GVat)^. 

Lemma 4.5. When F is a function on B, one has 

(4.29) \\{Vn{GVm?F)\\^ < C2\\{VnF)U + (^^^^^H^H- 

Proof. We use (4.25) and write 

(4.30) {Vn{GVn?F) = {VnAi) + (V^yls) + (FtvAs). 
With the help of Lemma 3.2 we find that 



(4.31) 



||(^iv^i)||oo ^<'^ -^\\VnG{VnF - {VnF))\1 
VlogiV 

(3-10) c 
- (logA^)3/2"^"°°' 



(3.18) r (3.8) 

(4.32) |KV^^2)||oo < -j==\\Vn[G{VmF)U^ < c\\{VnF)\U 
and that 

IK^^^3)||oo ^<'^ -=^\\Vjv{G{VnF) - [G{VnF)]o)\\oo 

VlogiV 

(4.33) 

(3.12) c 
- (logAr)3/2ll^ll~- 

Collecting (4.30)-(4.33), we obtain (4.29). □ 
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We now prove (4.12), (4.13). As a result of (4.24), (4.29) we see that for 
k>l and F a function on B, one has 

(4.29), (4.24) 

\\{VM{GVNrF)\\^ < C2\m{GVMf^'~'^F)\\^ 

(4.34) 

^ (logAf)3/2^1 " 

and by induction 

< 4\\ {VnF) IU + + C2ct' + ■■■ + 4-')l|i^l|oo 



<4\\{VnF)\\oo + 



(logAr)3/2H HOC- 

Keeping in mind that (Va^) = 0, we thus see that for k>l, 
|a.(2A. + l)|("y^-«)^^|(y.,(Gy.)-l)| 

(4.35) < alogN\A\\\{VN{GVNfh)\\^ 

(4-34) ac^lAI 

Together with (4.9), this proves (4.13) as well as (4.12) for odd n. When 
n = 2k, with A: > 1, we note that 

\aNm\ = ^^^\{VN,{GVN){GVNf'^''-'h)\ 

(4.36) =^}^\^GVN,VN{GVNf''''-'h)\ (by symmetry of G) 

(3.10), (4.24) 

<alog7V|A|||Gy;v||oo||V^(GV}v)'(^-i)l||^ < acK 

The proof of (4.12) is now complete. 

There remains to prove (4.14), that is, to analyze the large behavior 
of the even coefficients ai^{2k). To motivate the next lemma we recall that 

(4.37) a.(2) ^t|) ^iV,,GV,) ^\A\ Y^{V,GV,).. 

Lemma 4.6. There exists a function r(-) on N tending to at infinity 
such that 

(4.38) log N{VNGVN)z = \^^ + fN{z) for z e J 
with \fN{z)\<T{d{z,r)). 
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Moreover if one defines 

(4.39) TN = ^rj^ V g{{0,z)) so that limr^ = ^ by (1.9), 

one has the identity on B 

(4.40) {GVNflB - tn£{W)Ib = T-^MfN o vrz) + 

log A/ 

where for x = (y, z) G B, 

(4.41) |A;^|(a;)<_^log(iVMz,J^)). 

Proof. We begin with the proof of (4.38). We note that for z e J, 
X = {y, z) in B and x' = (y' , z') in B, one has 

log N{VnGVn). = 4t E Wiy)gix,x')W{y') 

' ' y,y'eA,z'eJ 

(4.42) = ^ t^(y)W^(y')E(<7(x'-x)-ff((0,z'-z))) 

^=°^ m E my)^(y')(fe7v(?/'-y,^)-^a(y'-y)) 

I I y,y'eA ^ ^ 

(4^5) \em 

where we have set 

/^(^) = m E "^{vW^y^Nky -y.z) iorz&j. 

The estimate in the second hne of (4.38) is now an immediate consequence 
of (1.11). This completes the proof of (4.38). 

We then turn to the proof of (4.40), (4.41). We write 

(4.43) {GVNfls = G{VnGVn) + G{VnGVn - {VnGVn)) = ai + aa. 
We know that 

(3.10) (3.10) 

(4.44) llaslloo < c\\VnGVn\\oo < c/logN, 
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and by (4.38) we find that 

1£{W) GIb 



ai 



+ 



1 



(4.45) 



2 |A| logiV logN 



TNSiW)lB + 



1 



loeiV 



where for x = {y,z) £ B we have set 



(4.46) 



"^^^^^-9 lAI lna/V"9W/V g[[0,Z-Z)). 



2 lAI logiV 2 1ogAf ■ 

' ' ^ \z'-z\<N 



We thus see that (recah Jy = {y} x J) 



|?-Ar(x)| < 

(4.47) + 



c 



V|A| 



+ 



logA^ 



1^*1 Jo - [Cljolol 



logiV 



E 9{{0,z'-z)) 



\z'-z\<N,z'(^.J 



(3.10), (3.12), (1.3) c 

Collecting (4.43)-(4.47) we have completed the proof of (4.40), (4.41). □ 

As a result of (4.37), (4.38) we see that 
(4.48) 



a 



a^(2) = ^£{W) + ^^f^iz)-, ^£iW). 

zeJ 



This proves (4.14) in the case k = l. To handle the case A; > 1, we will need 
to control the propagation of boundary effects corresponding to terms with 
Z-component close to the complement of J, when proving the convergence 
of aN{2k) = a-^{VN, {GVn)^''~^1) for N ^ oo. The next lemma wih be 
useful for this purpose. We first introduce some notation. 
We consider nonempty sub- intervals of J, 



Lo = d(/o, J \ /i) > 1, Li = dih,J \ h) > 1, 



(4.49) 
and define 
(4.50) 
as well as 

(4.51) Co = A X /o C Ci = A X Ii = C2 = A X I2 C B 

We have the following variation on Lemma 4.4. 
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Lemma 4.7. For F a function on B, one has 
\Ic,{GVn?F\\ 



oo 



(4.52) 

1,1, 1 1 A^ + 1) 



oo • 



Proof. By Lemma 4.4 we can assume that F vanishes on C2. We use 
the decomposition (4.25) of {GVnYF. We find that 

pco^illoo = \\IcoGVnG{VnF - {VnF))\\^ 

< \\lc,GVNlcMyNF - (F7vF))|U 

+ \\IcoGVnIb\c,G{VnF - {VnF))\\^ 



(4.53) 



(3.8), (3.9) 



< cViogiV||lc,G(FArF-(y^F))|U 

+ -=^log(^) \\G{VmF - {VnF))\U 

(3.11), (3.10) c C /iV+lV 

using the fact that F vanishes on C2 in the last step. In a similar fashion 
we find that 

||lco^2||oo = \\1coGVn[G{VnFM^ 

(3.10), (3.11) r 

< -^^\\1c,[G{VnFM^ 
+ -=^^||lB\Ci[G(y;vi^)]o||oo 

(3.9), (3.8) c /Af + 1\„ „ C „ „ 

< rrlog — — F U + 00, 

logiV \ Li J Lq 

where once again we have used that F vanishes on G2 in the last step. Finally 
we have 

||lCo^3||oo = \\IcoGVn{G{VnF) - [G(y7vi^)]o)||oo 
(3.8), (3.9) 



(4.54) 



^ ^ < c^/logN\\lcAG{VNF)-[G{VNF)]o) 

(4.55) 



(3.13), (3.12) c c /iV+ 1 

< -;-\\F\\oo + 77 log 



Li logN \ Lq 

using that F vanishes on G2 in the last step. 
Collecting (4.53)-(4.55), we obtain (4.52). □ 
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We now introduce the following sequence of possibly empty sub-intervals 
of J: 



Jfc = {i > 1; 1 + Ak[Ne-^^^] <i<N- 4/b[iVe-^/^]} 

(4.56) 

for A; > 0. 

There is some freedom in the above defin ition . The proof below would work 
with minor changes if one replaces ]\[e^^^°^^ by N^~^^ , with e^v — 0, and 
Eat log — > oo . 

Setting = A X Jfc, we find that for any k>l, when N > c{k), 
/ Sfc C Bk-i C • • • C 5i C S 
and that with the notation (4.39) 

\\lB,iiGVNriB-irNSiW))'lB)\\^ 
ifc-1 

m=0 

(4.57) 



B 



I oo 



jfc-1 



<Y,{rN£{W)?-^'^^'^UB^AGy?"'{{GVN?lB-rN£{W)lB)\\^. 

m=0 

We set Fat = {GVNflB-TN£{W)lB, and now want to bound (GVat)^™ x 

-^Af||oo with the help of (4.52), when < m < fc. To this end we introduce 
Jm with a similar definition as in (4.56), simply replacing Ak by 2m + 2, so 
that Jrn+i Jm Jm play the role of /q ^ ^ h in (4.49). We note that 
d{Jm+i, J \ Jm) and d{Jm, J \ Jm) are bigger than cNe~^^°^^ . As a result 
the expression inside the parenthesis after C5 in (4.52) is smaller than 



log A/ \ 



(N + l 



\2 



Af2 ) - ^/b^' 

It thus follows that for < m < A; 

(4.52), (4.24) 

, , \\Ib^AGVn)^^Fm\\oo < c41b^{GVm)^^^-'^F^\\^ 
(4.58) 

^m-l 

+ ^=iC6||FAr||oo 

VlogiV 

and by induction 

< ctwiB.F^Woo + ^=^(cr ' + C4cr-' + • • • + cr')ii^^iic 

Vlog A/ 
VlogiV 
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Coming back to the last line of (4.57), we see that for A; > 1, > c{k), each 
term under the sum, thanks to the above bound and (4.39), is smaller than 



c(A:)(||1b,F^||oo + 



1 



(4.59) \\lB,iiGVNf^lB - {rN£iW))'lB)\\^ < c(fc) ( ||1biF7v||oo + 



||F/v||oo)- Hence we see that for /c > 1 and N > c{k), 



1 



where we have used the bound ||-F/v||oo ^ c, which follows from (4.24) and 
(4.39). Note that the form of the correction term (log A^)^^/^ in (4.59) mainly 
reflects our choice for the intervals Jk in (4.56). In view of Lemma 4.6, we 
also find that 



\Ib,Fn\ 



(4.40) 

< 



(4.60) 



1b 



G{fN ° TTz) 



logA^ 

+ lllfii^lloo 



oo 





as N ^ oo. 



where we have used the bounds in the second line of (4.38) and (4.41) to 
conclude in the last step. Since r^v converges to ^ [see (4.39)], we can infer 
from (4.59), (4.60) that 



A 



dcf 



(4.61) 



2tt 



N 







for each k> 1. 



We can now use this estimate to study the asymptotic behavior of aAr(2( 
1)) as goes to infinity. Indeed one has 



aN{2ik + !))-( ^£{W) ] aNi2) 



(4.62) 



a- 



logN 



VN,GVN(^iGVN?''lB- (^-^£{W)^ 1b)) 



N 

= h + h, 

where in the last step, using the symmetry of G, we have set 



Ji = Q 



logA^ 



: a- 



N 
logN 



GVN,VNlBA{GVNf''lB 



2tt 



SiW)] 1b 



N 



GVn, Vn1b\b, ( {GVNf^lB - ( ^£{W) ) 1b 



27r 



We then observe that 



(3-10) , (4 61) 

/i <ca log 7V||GViv Hoc II II ooA,.,,v < c'aAk,N -^ 
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and that 

(4-24) \B\Bk\ 

h < alogiV||GV^||oo||^7v||oo xc(fc) X ' ^ " 



N 



We have thus shown that 

(4.63) hm aAr(2(A; + l)) - f^£:(T^)^ aAr(2) =0 forany/c>l. 
TV \ zvr / 



Combined with (4.48) this completes the proof of (4.14) and hence of The- 
orem 4.1. □ 

Our next objective is to study the convergence in distribution of the ran- 
dom fields ( j^g'j^ )^gg2 and {j^^)yez^ as N goes to infinity, where we recall 
the notation from (4.1), (4.2). The task is simplified by the fact that we have 
already proved Theorem 4.2: we only need to investigate the convergence in 
distribution of these random fields at the origin. We now focus on the case 
where A = {0}, and Jq plays the role of B. We further define 

(4.64) Vn{x) = -^1jo{x), V^{x) = —Vn{x) for x G 

log N 

and set 

(4.65) c^ = -A_Co,N, C'^ = ^£'o,N. 
Just as in (4.7), (4.10), we know by Theorem 2.1 that 

(4.66) E[exp{z£Ar}] =exp< ^^aAr(n)z" > for \z\ < Fat in C 

with r^r > and where we have set 

(4.67) aN{n) = UN{VN,{GVNT~^l) forn > 1 
and that 

(4.68) E[exp{z£'j^}] = expi ^S'^(n)z" i for \z\ < in C 
with r'j^ > and 

(4.69) a'M = u'AVl,, {GVl,r-H) = l^l^-^a^{n) for n > 1. 

a 

The heart of the matter for the proof of Theorem 4.9 below lies in the control 
of the large N behavior of the sequence ai\i{n), n > 1. 
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Theorem 4.8. 

(4.70) 0<aN{n)<ac^ forn>l,N>l, 

(4.71) for any n>\ limaAr(n) = a( — 

N yvr 

We first explain how Theorems 4.2 and 4.8 enable us to infer the con- 
vergence in law of the appropriately scaled fields Cy^N, y G Z^, and >C^^, 

Theorem 4.9. 

As N goes to infinity, {\^^)y^7? converges in distribution to a flat 

(4.72) random field with constant valued distributed as with R as in 
(4.18). 

c 

As N goes to infinity, { ^^f )y^'p converges in distribution to a flat 



(4.73) 



random field with value 1. 



Proof (assuming Theorem 4.8). A repetition of the arguments used 

in the proof of Theorem 4.2 shows that Cn converges in distribution to a 
nonnegative random variable C with Laplace transform 

(4.74) ^[exp{-A£}] = expl 1 for A > 0, 



l + SA/vr, 

so that by (1.30), C is BESQ'^(a, ^)-distributed, that is, has the same dis- 
tribution as V"^ in the notation of (4.18). 

Moreover we know from Theorem 4.2 that for any y G Z^, when N goes 
to infinity, jy i^y.N — ^o,n) converges to zero in distribution, and (4.72) 
follows. 

In the case of (4.73) we note instead that the arguments used in the proof 
of Theorem 4.2 now show that Cj^ converges in distribution to a nonnegative 
random variable with Laplace transform e~^, A > 0, that is, to the constant 
1. Since by Theorem 4.2, j^{^'y n n) converges in distribution to zero 

for any y G 1?, we obtain (4.73). 

Proof of Theorem 4.8. We now write || • ||oo for the supremum norm 
on B = Jq, and the ^inear operators we consider are restricted to functions 
that vanish outside B. The fact that ON^n) is nonnegative is plain; see (4.64), 
(4.67). Moreover the right-hand inequality in (4.70) is a direct consequence 
of (3.8). This proves (4.70). 

We now turn to the proof of (4.71). For A; > 1, we introduce B^ = {0} x 
Jfc, with Jfc as in (4.56), so that for any A; > 1, and N > c{k), (j) ^ B^ $1 
B^^i • • • C 1?]^ C S. In a much simpler fashion than (4.59) we now find 
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that [see (4.39) for notation] 

(4.75) \\l^^{{GVMfl^ - (2r^)h5)|L < c(/c) (jjl^ F^|U + ^j^) , 

where we have set 

Fn = (G'Vjv)lg - 2tn1§ = GVn - 2tn1^. 
It already follows from the definitfons of V/v and tn in (4.64), (4.39) that 

(4.76) lim||l^^F7v||oo = 0. 
We thus see that for /c > 1 , 

limS;v(A:) = lima^j^{VN, (GVn)''-'!) 

(4.77) 

/ /3' 
= a\^2 lim tnj = al — 



and this proves (4.71). □ 

Remark 4.10. (1) There is an important connection between random 
interlacements at level u and the structure left by a random walk on a large 
torus {Z/NZ)'^ (here d = 3), at a microscopic scale of order 1 (see [23]) or 
even at a mesoscopic scale of order N^~^, with < e < 1 (see [22]) when 
the walk runs for times of order uN"^. This naturally raises the question 
whether the above limiting results might also be relevant for the field of 
occupation times left close to the origin, by continuous time simple random 
walk with uniform starting point on a large two-dimensional torus (Z/A^Z)^, 
at times of order aA^^logA^ {=uj\fN^) or at much larger times u'j^N^. Let 
us incidentally point out that the time scale aiV^ log N is much smaller than 
the cover time of the torus which has order ^A^^(log A^)^; see [4]. 

(2) We can consider the discrete occupation times £x,u7 x G Z'^, u>0 [see 
Remark 2.4(5)] and define £y^N and ^'y^^, for y e 1? , N > 1, as in (4.2), 
simply replacing La;,^, by 4,„. 

£' -£' 

in distribution to a centered Gaussian field which vanishes at the origin. 
However, this limiting field is different from (V'?/)j/ez2 in (4.20). 



Theorems 4.2 and 4.9 enable us to show that ( ^''J- — ^^ )y^-^2 converges 



The heart of the matter lies in the fact that for any W{-) as in (4.3), when 
N tends to infinity, 

y^o,cy2 VF(v)^== converges in distribution to a Gaussian variable 
(4.78) ^^^^ 

with variance £{W) — ^y^i2 W{y)'^ 
[in the case of C'y ^ the limiting variance instead equals i5(M^)]. 



36 



A.-S. SZNITMAN 



Indeed it follows from (2.27) that for real z and > 1, 



E 



(4.79) 



= E exp J^(l-e-^^(j')/^^)4,^ 
and for |z| < r and N >c, we can use Taylor's expansion and write 



1 -e" 



N 



Nw 



^(l + e,(z,iV)), 



N 



\z\ < r. 



where \ey{z,N)\ < ^ and Huin ey{z, N) = 0, for each y G 

Inserting the above identity in (4.79) shows that for \z\ < r and N > c, 
the left-hand side of (4.79) equals 



(4.80) E 



By the end of the proof of Theorem 4.2 we know that for \z\ < c, 



limE 

N 



exp<^ z ^ 



W{y) 



''N 



:exp<j -f(M^) 



A straightforward uniform integrability argument combined with (4.73) and 
(4.80) shows that for real z with j^l < c, 

W{y) ^, .2 



(4.81) limE 

N 



exp< z 



E 



ye2 



: exp 



Similar arguments as in the proof of Theorem 4.2 now yield (4.78). 

Note incidentally that Theorems 4.2 and 4.9 are not quite sufficient to 



study the limit in law of ( 



. As shown by the above proof 



[see in particular (4.80)], to handle this case we would in essence need a 
limiting result for the joint law of the two random fields that appear in (4.17) 
and (4.72). □ 



5. Occupation times at high level u. In this section we relate occupation 
times at a high level u of the random interlacements with the d-dimensional 
Gaussian free field. The limit u — )• oo bypasses the obstructions present when 
one considers a fixed level u\ see Remark 2.4(2). Our main result appears in 
Theorem 5.1. It has a similar flavor to (4.20) of Theorem 4.2 and (4.73) of 
Theorem 4.9. Moreover it can rather straightforwardly be extended to the 
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case of random interlacements on transient weighted graphs; see Remark 5.2. 
However, we keep the set-up of Z'^, d>3, for the main body of this section, 
not to overburden notation. 

We consider on an auxiliary probability space 

. ^x,x G Tj'^, a centered Gaussian field with covariance function 

^ • ^ E[ip.^ip^,]= g{x' - x) + g{0) - g{x) - g{x'), for x,x' GZ"^. 

This field has the same distribution as the field (73; — 7o)a;eZ'* of increments 
at the origin of the d-dimensional Gaussian free field, (7a;)xeZ''; ^^^^ is, 
the centered Gaussian field with covariance function E^y^Jx'] = 9{x,x'), for 

We can now state the main result of this section. 



Theorem 5.1. As 00 



2u )a:eZ'* Converges in distribution to the Gaussian random 



(5.2) ^ 

field {iPx)xGZd ^ri (5.1) 



and 
(5.3) 



—Lx^u I converges in distribution to the constant field equal 
" ' J xezd 
to 1. 



Proof. We follow the same strategy as in the previous section, the 
problem is, however, much simpler now. It clearly suffices to prove (5.2) 
and (5.3) with u replaced by a sequence un such that 

(5.4) UN > 1 for > 1 and limuN = 00. 

We thus consider a function V on Z"^ such that 

(5.5) V is finitely supported and V{x) = 0. 

We define 

(5.6) LN=y] -j^V{x)Lx,uj, for N>1. 

x&.<^ ^ 

It follows from Theorem 2.1 and Remark 2.2 that for some fixed r > 0, 

(5.7) E[exp{2;LAr}] =exp< ^^CAr(n)2;" > for \z\ < r in C, 

where we have set for n,N > 1, 

(5.8) CN{n) = 2-^/\'-"/\v,{GVr-H). 
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In view of (5.5) we have 

(5.9) C7v(l) = 0. 
By (5.4), (5.8) it is also plain that for n > 2, 

(5.10) \cN[n)\<c{yr 

with c(y) a positive constant depending on V and d, by our convention. 
Moreover we find that 

(5.11) li^CAr(?i)=0 for n > 2 



and 



(5.12) CM{2) = \{V,Gvf-'^d'-'^\E 



The same arguments as in the proof of Theorem 4.2 show that 



L jv converges in distribution to a centered Gaussian variable with 
variance E[{Y^^^V{x)ipx)'^]- 



Since V in (5.5) and um in (5.4) are arbitrary, claim (5.2) follows. We then 
turn to the proof of (5.3). It follows by (2.21) that for A > 0, 



E 

(5.14) 



exp<' - —Lq uj. 

UN 



exp 



A 



l + g{0){\/uN) 
as — 7- oo. 

This shows that ■:^Lq^ui^ converges in distribution to the constant 1 as 
goes to infinity. Since due to (5.2), for any x G l/, :;^{Lx,uiq — -^o,ujv) tends 
to zero in distribution, as goes to infinity, our claim follows. □ 

Remark 5.2. (1) The results of the present section can straightfor- 
wardly be extended to the set-up of continuous time random interlacements 
on a transient weighted graph E, as we now explain. We keep the same 
notation and assumptions as in Remark 2.4(4). We introduce a base point 
xq £ E. In place of (5.1) we consider 

fx,x £ E, a centered Gaussian field with covariance function 
E[(px(Px>]=g{x,x')+g{xo,xo)-g{xo,x)-g{xQ,x'),x,x' £E, 

where g{-,-) now stands for the Green density. 

This field has the same distribution as the field of increments {jx — 
lxo)xeE, of the Gaussian free field (7x)xe-B attached to the transient weighted 
graph, that is, the centered Gaussian field with covariance function E[yx^x'] = 
g{x, x'), X, x' S E. 
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With the help of (2.25), (2.26), the arguments employed in the proof of 
Theorem 5.1 now show that as n — t- oo. 



(5 16) V^M )^S-E converges in distribution to the Gaussian random 

field (¥?x)a;eB 

and that 

(5.17) [ —Lx.u\ converges in distribution to the constant field equal to 1. 



^ / X6-B 



(2) In the case of the discrete occupation times £x,u^ ^ G ^'^i u>0, the 
same arguments used in Remark 4.10(2) show that when u — t- oo, 

(5.18) f _x^u_j£fi\ converges in distribution to (^'x)a;6Z'*5 
V v2u 



where {t^x)xeZ'^ centered Gaussian field vanishing at the origin such 

that for any V as in (5.5), one has in the notation of (5.1) 



(5.19) E 



Moreover looking at the Laplace functional, one sees with the help of (2.27) 
and (5.3) that for u— )• oo, 

(5.20) 

~(-x,u I converges in distribution to the constant field equal to 1. 
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